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Comparative analyses of protein-protein interaction networks play important roles in the under¬ 
standing of biological processes. However, the growing enormity of available data on the networks 
becomes a computational challenge for the conventional alignment algorithms. Quantum algorithms 
generally provide greater efficiency over their classical counterparts in solving various problems. One 
of such algorithms is the quantum phase estimation algorithm which generates the principal eigen¬ 
vector of a stochastic matrix with probability one. 

Using the quantum phase estimation algorithm, we introduce a quantum computing approach for 
the alignment of protein-protein interaction networks by following the classical algorithm IsoRank 
which uses the principal eigenvector of the stochastic matrix representing the Kronecker product 
of the normalized adjacency matrices of networks for the pairwise alignment. We also present a 
greedy quantum measurement scheme to efficiently procure the alignment from the output state of 
the phase estimation algorithm where the eigenvector is encoded as the amplitudes of this state. 

The complexity of the quantum approach outperforms the classical running time. 


I. INTRODUCTION 

Comparative analyses of protein-protein interaction 
(PPI) networks play important roles in the understand¬ 
ing of biological processes. Alignments of PPI networks 
drawn from different species provide invaluable infor¬ 
mation to catalog conserved network regions and iden¬ 
tify functional similarities across species. Using different 
formulations, many network alignment algorithms have 
been proposed such as the ones in 0-9 and others @. 
However, the growing enormity of available data on the 
networks raises computational challenges for implemen¬ 
tations of these algorithms. In analogous to ranking al¬ 
gorithms, e.g. PageRank [7|, one of the global alignment 
algorithms, IsoRank [§}, uses the intuition that the score 
of aligning two nodes should depend on the alignment 
of their neighbours in the global alignment of two PPI 
networks. It formulates the network alignment as an 
eigenvector problem where coefficients of the principal 
eigenvector of a stochastic matrix, the stationary state, 
represent the functional similarity scores between pairs of 
the nodes. Using a greedy algorithm, it then generates 
the network alignment from the eigenvector. In the case 
of multiple networks, IsoRank is applied to every pair 
of the networks and a global alignment is retrieved from 
the pairwise network alignment results [9j. While Iso¬ 
Rank is one of the successful algorithms and fast enough 
to handle the alignment of large sparse graphs 0 , the 
exponential scaling of its running time with the number 
of networks still impedes applications of the algorithm to 
multiple networks (the running time of the algorithm is 
0(E m ) where E is the number of edges in a network and 
m is the number of networks) 0 - 

In 0|) it has been discussed that aligning networks 
on quantum computers may provide greater efficiency. 
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In particular, it is shown that the quantum phase esti¬ 
mation algorithm can be used for stochastic matrices to 
find their principal eigenvector with the success probabil¬ 
ity one. However, it has not been shown how to procure 
a solution from the quantum state representing the prin¬ 
cipal eigenvector. It is also not shown how to simulate 
ranking matrices which are not Hermitian in most of the 
cases in the phase estimation algorithm. 

Here, we first show explicitly how the eigenvector of a 
stochastic matrix can be obtained by using the quantum 
phase estimation algorithm by following mostly Ref. [0. 
To construct the alignment from the final quantum state, 
we then give a greedy quantum algorithm based on a 
quantum measurement scheme. In addition, we present 
an approximate simulation approach for non Hermitian 
ranking matrices and present simple examples with nu¬ 
merical results which can be experimented on quantum 
computers (The capacities of the current quantum com¬ 
puters are very limited. Thus, they can only run for small 
sized problems.). In the end, the complexity analysis 
shows that the multiple network alignment by this quan¬ 
tum approach requires o(m x poly{log{ |U|))^ computa¬ 
tional running time for sparse matrices while 0(m |U| 2 ) 
for dense matrices. This is an exponential speedup over 
the classical running time. 


II. CLASSICAL ISORANK 

A PPI network is generally represented as an undi¬ 
rected graph, G(V,E) where the set of nodes V repre¬ 
sents the set of proteins and the set of edges (vt, Vj) € E 
describes interactions between proteins Vi and Vj . To ob¬ 
serve conserved similarities across species, PPI networks 
for different species are comparatively analyzed by maxi¬ 
mizing an objective function to indicate correspondences 
between the nodes. One of the most common strate¬ 
gies to observe the similarities between networks is the 
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network alignment. Network alignment is the process to 
globally compare two or more networks to identify the re¬ 
gions where these networks are similar and dissimilar. In 
general terms, network alignment can also be considered 
as a more general variant of the subgraph isomorphism 
problem which is an NP-complete problem of determin¬ 
ing whether a finite larger graph contains a smaller graph 
as an exact subgraph [13] . However, in the network align¬ 
ment, it is important to find the most similar and dissimi¬ 
lar parts between two graphs even if the larger graph does 
not contain the smaller one as an exact subgraph. jufl] 
Singh et al.Q have presented a global alignment algo¬ 
rithm, viz. IsoRank, by using the intuition that the score 
of aligning two nodes should depend on the alignment of 
their neighbors. This intuition is formulated as follows: 

Rij= ^ T, \n(u)\\ N (v)\ R uv ’ (1) 

u£N(i)veN{j) 1 v n 

where, N(a) is the set of neighbors of the node a; |iV(a)| 
is the size of this set; V± and Vi are the set of the nodes 
for the networks G\ and G 2 ; and i £ V i and j £ Vi. 
R defines the functional similarity matrix whose station¬ 
ary state is used to find the solution for the alignment 
problem. Eq.Jl]) can be rewritten also in matrix form as: 

R = AR, 

au j]W v] = (PM if (*>“) e and t7» e E 2. 

’ ’ l 0 otherwise. 

( 2 ) 

In the above equation, A is a |Vi11Va|x|Vi11V 2 I stochastic 
matrix and can be defined from the Kronecker product 
of the column-wise normalized adjacency matrices of the 
input graphs: A = A\ ® A 2 , where Ai is the normalized 
adjacency matrix for the graph Gi and also a stochastic 
matrix. Eq.© describes an eigenvalue problem where 
the principal eigenvector corresponding to the eigenvalue 
one of the matrix A is the stationary distribution of the 
random walk on the Kronecker product graph. On clas¬ 
sical computers, this equation can be solved by using 
different iterative methods such as the power iteration. 


III. QUANTUM PHASE ESTIMATION 
ALGORITHM FOR NETWORK ALIGNMENT 

In this section, we briefly explain the main intuition 
of the quantum phase estimation algorithm (for an un¬ 
familiar reader, we recommend Ref. [151]) and then show 
how the stochastic matrices can be used in the phase es¬ 
timation algorithm by following mostly Ref.yjl]. Then, 
we discuss how to simulate a non-Hermitian stochastic 
matrix on quantum computers by using the closest Her- 
mitian matrix. In the final subsection, we show how to 
incorporate some other information used in the PPI net¬ 
works. 


A. Quantum Phase Estimation Algorithm 

For a given approximate eigenvector encoded as the 
amplitudes of the quantum state | Hj) and the eigenvalue 
equation U |/u}_= \g,j)] the phase estimation algo¬ 

rithm (PEA) [1|| tries to find the phase <f>j in this equa¬ 
tion. PEA mainly requires two quantum registers, \regi) 
and \regi), consisting of sufficient number of qubits to 
hold the eigenvector and the phase, respectively. In the 
initial setting, \regl) is set to zero state and |re<? 2 ) is 
assigned to hold a vector which is the best known ap¬ 
proximation of With the help of quantum Fourier 
transform and the sequential controlled unitary opera¬ 
tions, \reg\) becomes holding the Fourier transform of 
the phase. Then, the application of the inverse quantum 
Fourier transform turns \regl) into the binary value of 
the phase: \regl) = | <f>j). Consequently, the value of the 
phase is obtained by measuring \regl) in the standard 
basis. Here, if the unitary operator U is the time evolu¬ 
tion operator of a Hermitian matrix H 1 U = e t2TH , then 
one also obtains the eigenvalue of H. 


B. Application to Stochastic Matrices 

The success of the phase estimation algorithm is di¬ 
rectly related to the closeness of the input vector to the 
actual eigenvector. This can be defined by the dot prod¬ 
uct. The dot product of an equal superposition state and 
a vector is the normalized sum of the vector elements. On 
the other hand, the eigenvectors of a stochastic matrix 
has the property that the sum of the vector elements is 
one for the principal eigenvector and zero for the rest of 
the eigenvectors. In [12[, it has been showed that when 
PEA is given an equal superposition input state, it then 
finds the principal eigenvector of U and so the principal 
eigenvector of the stochastic matrix H with the success 
probability equal to one. 

In our case, we find the principal eigenvector of the ma¬ 
trix A which is the Kronecker product of the column-wise 
normalized adjacency matrices (normalized in the sense 
to have stochastic matrices) for the input networks: i.e., 
A = A\ <g> Ai 0 ■ ■ ■ <g) A m . Because of the Kronecker 
product, A defines a separable system and so quantum 
circuits for each Ai can be constructed separately. This 
eases the difficulty of finding quantum circuits for the 
simulation. However, in general, quantum computing is 
based on unitary gates associated with time evolution op¬ 
erators of Hermitian quantum systems. This dictates the 
stochastic matrix used in the algorithm to be Hermitian, 
in which case the principal eigenvector is already known 
to be a vector of all ones. In the following subsection, we 
shall describe how to approach non-Hermitian matrices. 
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C. Simulation of Non-Hermitian Operators 


A matrix A is called positive if the matrix elements 
Aij > 0 and non-negative if A^ > 0. It is normal if A* A— 
AA 1 = 0, where A t describes the conjugate transpose of 
A. Any matrix A can be decomposed into a Hermitian 
and a skew-Hermitian matrices as: 

A = H + S = ±(A + A) + 1(A-A% (3) 

where H = 4(A + A*) and S = ^(A — A) define the 
nearest Hermitian and skew-Hermitian matrices to A, re¬ 
spectively E3- The eigenvalues of H are all real and the 
eigenvalues of S have only imaginary parts. Moreover, 
when A is a normal matrix, there are a few additional 
useful properties: 

• H and S commute: [H, 5] = HS — HS = 0. 

• Since A A 1 = A^A, H and S have the same eigen¬ 
vectors. 


• The imaginary part of the eigenvalues of A are 
equal to the eigenvalues of S, and the real parts 
are equal to the eigenvalues of H. 

Because of the last property, one can simulate normal ma¬ 
trices and their corresponding non-Hermitian operators 
on quantum computers by using two separate registers to 
obtain the imaginary and the real parts of the eigenvalue 
individually. In that case, one uses two unitary opera¬ 
tors U\ = e lH and U 2 = e S for the simulation. (Note 
that U 2 is a unitary matrix because the exponential of 
a skew symmetric matrix is a unitary matrix.). How¬ 
ever, if a stochastic matrix is normal, it turns out that 
it is also doubly stochastic: i.e., its left and right princi¬ 
pal eigenvectors are known to be a vector of all ones with 
the eigenvalue one. Therefore, instead of an approximate 
normal matrix, we shall use the closest Hermitian matrix 
H = i(A +At) in our simulations. H defines a non¬ 
negative matrix. If it is also irreducible, then there is a 
well known theorem stating that H has a dominant posi¬ 
tive eigenvalue equal to the spectral radius, p(H ), of ma¬ 
trix H and corresponds to a real nonnegative eigenvector 
unique up to a scalar multiple [l^|. The value of p(H) is 
bounded as: mirij hij < p(H) < maxj hij , where 
hjj are the matrix elements of H. Using the stochastic 
property of A, we also find p{H) > 0.5. PEA outputs 
jth eigenvector with the probability: 


1 

Vn 


J2 ('I IAj) 


( 4 ) 


This is the normalized sum of the coefficients of the j th 
eigenvector and is equal to one if | /x) ■ is the principal 
eigenvector and zero for the rest of the eigenvectors. 
Therefore, the success probability of PEA for stochastic 
matrices is one. Moreover, since H is the best Hermi¬ 
tian approximation to A, the principal eigenpair of H 


close to the ones of A. Hence, we expect the sums of the 
coefficients of the eigenvectors of H to be very similar 
to the ones of the eigenvectors of A and so the success 
probability to be very close to one. 


D. Incorporation of Other Information 


As done in IsoRank [g], one can include further infor¬ 
mation, e.g. BLAST scores, into the quantum model as 
well in the following form: 


A = H + B, 


( 5 ) 


where B and H are to be assumed to commute: BH — 
HB = 0. Therefore, the time evolution can be written 
as: 


e iA = e iH+B = e iH e iB 


( 6 ) 


Note that the above equation does not change the spar¬ 
sity of the matrix A; hence, it is still sparse and the 
evolution operator and the corresponding quantum cir¬ 
cuit require polynomial time for the implementation fl9| 
(see Sec IVIl for the complexity analysis). 
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E. Summary of the Method 

• Prepare to quantum registers as: |re< 7 i )=0 
and \reg 2 ) = -j= Yli |i) defining an equal 
superposition state. 

• Find H = A+ 2 Af as the closest Hermitian 
approximation to the stochastic matrix A. 

• Find the quantum circuit simulating H by 
using one of the methods described in [19!, 
[20| . The number of gates in the resulted 
circuit is polynomially bounded w.r.t. the 
number of qubits for the sparse matrices. 

• Apply the phase estimation algorithm with 
the above initial settings. 

• For stochastic matrices; the PEA results the 
correct answers with probability one: the 
probability of seeing jth eigenvector is given 
by the normalized square of the sum of its 
coefficients: 


1 

Vn 


N 




(7) 


Since the sum of the coefficients are one for 
the principal eigenvector and zero for the 
rest, PEA results the principal eigenvector 
with probability one. 

• For the closest Hermitian matrix U; the 
probability is not one. However, it is still 
expected to be very close to one since H is 
a nonnegative matrix with a principal eigen¬ 
value and a positive eigenvector (this is due 
to the Perron-Frobenious theorem) close to 
the eigenvector of the original stochastic 
matrix. 


we cannot apply classical matching algorithms directly. 
Hence, engineering the order of the quantum registers in 
the measurement, we describe the following greedy strat¬ 
egy for the alignment: Consider the eigenvector consists 
of two registers as |x) |y), where each register represents 
a network: 

1. Apply a conditional measurement: when the first 
register is Xi, measure the second register. For each 
measured Xi, this generates a yj; therefore, XiS are 
matched to yjS. Note that the measurement out¬ 
come is determined by the conditional probability 
(The probability is conditional when additional in¬ 
formation such as BLAST scores is included. Oth¬ 
erwise, it is disjoint.): the probability of measuring 
xjj in the second register while the first register is 

Xi. 

2. If there are duplicities and still unmatched nodes; 
then apply a second type of conditional measure¬ 
ment: if the second register yj , measure the first 
register. This time for each yj an Xi is obtained. 

This step is only useful when additional informa¬ 
tion is incorporated (A + B is used) because the 
system is otherwise separable as A = A\ ® A 2 and 
the conditional measurement outcome will be the 
same as the outcome obtained in the first step. 

3. We combine these two different measurement out¬ 
comes obtained in the first and the second steps 
and begin matching from nodes whose scores are 
the highest. 

4. If there are still unmatched nodes, statistical infor¬ 
mation about the other possibilities obtained dur¬ 
ing the measurement is used to match the remain¬ 
ing nodes. 

Since the main intuition of the algorithm gives a higher 
similarity score to the nodes whose neighbours have 
high scores; in the matching, the neighbours of the first 
matched nodes are given priority so as to generate a so¬ 
lution which is also connected. In the end, we choose the 
largest connected component as the best solution to the 
alignment. 


B. Matching for Multiple Networks 


IV. EXTRACTING NODE MAPPINGS FROM A 
QUANTUM STATE 

A. Matching for a Pair of Networks 

Generating a discrete solution from the final quantum 
state is known to be just solving a maximum weight 
matching problem. However, since fully obtaining a 
quantum state requires exponential time complexity [l5[ , 


In the case of multiple networks, as done in Q an ob¬ 
vious approach is to align each possible pair of networks 
and derive a common solution for the multiple networks 
from these pairwise alignments. However, the complex¬ 
ity of this approach grows exponentially with the number 
of networks. Here, we shall follow a different approach 
which can be generalized easily to any number of net¬ 
works: Assume we have three networks Gi, G 2 , and G 3 
with nodes {ai... a mi }, {61 ... b m2 } and {ci... c m3 }, re¬ 
spectively. Let also \regcQ, | feffG 2 }, and |fe<?G 3 ) repre¬ 
sent the networks Gi, G 2 , and G 3 , respectively: 
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V EXAMPLES 


• Measuring |re<?Gi) alone and \regG 2 ) and \regc 3 ) 
together, we draw the conditional probabilities to 
see one of the nodes, Gtj, in |re< 7 Gi) and bjCk in 
\regG 1 )\f’egG 3 ) ■ In other words, when |re<7Gi) is 
ai, the probability to see bjCk in ^egG^regcs) is 
obtained. This results in a matching of a; — bjC k - 
Note that in the real implementation on a quantum 
computer, one just assigns the nodes as a result of 
the measurement outcome. 


— We match nodes initiating from the largest 
probability, or the most commonly seen mea¬ 
surement outcome. 

— As done in the case of two networks, the prior¬ 
ity given to the neighbours of the first matched 
nodes. 


• If the measurement is no longer provide valuable in¬ 
formation to match further nodes, then we measure 
| regc 2 ) alone and |re< 7 Gi) an d | regG 3 ) together. 
This gives the probabilities to see one of the nodes, 
bj, in | regc 2 ) and a,iC k in \reg Gl )\reg G3 )- Then we 
combine this measurement result with the previ¬ 
ous measurement result to match the unmatched 
nodes. 

One can also go further and draw probabilities for 
Cfc — a,ibj and combine them with the previous results. 
However, while different measurement settings increase 
the statistical confidence in the measurement results, it 
increases the complexity of the algorithm. As noted be¬ 
fore, if \reg 2 ) is separable at the beginning: i.e., |reg 2 ) = 
|regGi) ® I re 9G 2 ) ® ® \ re 9G m )i then the other mea¬ 

surements also produce the same result as the first mea¬ 
surement setting, Gu — G 2 j ■ ■ ■ G m k . Thus, only the first 
measurement setting is used to conclude the matching of 
the nodes. 


V. EXAMPLES 


A. Example-1: Alignment of Network Pairs 


In the case of the alignment of two networks G 1 -G 2 , 
the possible probability outcomes are shown in Fig 2(c)} 
where it is first assumed that \regc 1 ) = ai, then prob¬ 
abilities in the collapsed state are found for each bj. 
Note that since the collapsed state is not normalized, 
it represents conditional probabilities: e.g., the proba¬ 
bility of |re<7G 2 ) = b 2 when \regG 3 ) = ai- The same 
approach is also applied to G 2 -G 3 and the probabilities 
are shown in Fig 2(a)| The matching algorithm applied 
to the outcomes in Fig |2(a)| and Fig 2(c) procures the ex¬ 
act alignments for the both pairs of the networks. To 
make the method understood well, we shall go through 
the alignment for G\ and G 2 by assuming PEA is run 
statistically enough times with the same setting: The 
probabilities for |re<?Gi) in the final state of PEA are 
[0.1176,0.1176,0.2647,0.1176,..., 0.1176,0.0294] where 
we see as with the highest probability. When 
\regGi) = 03 , the quantum state collapses to the follow¬ 
ing unnormalized state : [0.1213,0.3638,0.2425,0.2425] 
which produces to the following normalized probabilities 
[0.0556,0.5000,0.2222,0.2222]. From this state we see 
\regG 2 ) = b 2 with the highest probability. Hence, 03 
is matched to b 2 . In the next possible step, |re< 7 Gi) 
can take the same probabilities for different states. In 
the measurement, we give the priority to the neigh¬ 
bours of the matched nodes. When |re< 7 Gi) = ai, which 
is a neighbour of 03 , the quantum state collapses to 
[0.08080.2425,0.1617,0.1617] producing the normalized 
probabilities: [0.0556,0.5000,0.2222,0.2222], Since b 2 is 
already matched, the second highest possible outcome is 
either 63 or 64 . Therefore, ai is matched to either 63 or 
64 . In the final round, we match a 2 to the remaining one 
of 63 and 64 . The same but unnornralized probabilities 
are also shown in Fig|2(c) 


Fig|2(b)| and Fig 2(d) also show the probability out¬ 


comes for the alignment of the same networks: G\-G 2 
and G 2 -G 3 ; however, in these figures, instead of A, ap¬ 
proximate Hermitian matrices found by H = 1/2(A+A.1) 
are used. While the probabilities slightly differ in com¬ 
parison to Fig |2(a)| and Fig 2(c) they still conclude with 
the same alignments. 


Because of the computational difficulty in simulating 
quantum computers on classical computers (The com¬ 
plexity grows exponentially with the number of qubits 
involved in the simulation. PEA requires also sufficient 
number of qubits on the first register to hold eigenvalue.), 
here we only aim to show the capability of the algorithm 
when a fully functional quantum computer is available. 
Therefore, we shall use merely trivial examples given in 
Fig H] for which one can see optimal solutions easily. For 
the large non-trivial networks, since both IsoRank and 
the quantum approach are based on the principal eigen¬ 
vector; we expect the success of the algorithm to be sim¬ 
ilar to IsoRank even though the matching algorithm de¬ 
fined here is different than the one in IsoRank. 


B. Example-2: Alignment of Multiple Networks 


As an example for the multiple network alignment, we 
use Gi, G 2 and G 3 shown in FigQ] and follow a simi¬ 
lar method to the alignment of two networks: First, we 
assume |re< 7 Gi) = a*, then find the probability of see¬ 
ing \regG 2 regG 2 ) = bjC k in the unnormalized collapsed 
state which represents the conditional probabilities. The 
probabilities are shown in Fig 3(a) where the exact ma¬ 
trix is used. Fig 3(a) represents the probabilities when 
the Hermitian matrix H = 1/2(A + A') is used. Giving 
the priority to the neighbors of the first matched nodes 
as explained in Sec II V Bl the exact same matching can 
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(a)Gi 


(b)G 2 


(c)G 3 


FIG. 1. Example Networks 




(a)Probabilities for matching of the nodes of G2 and G3 (b)Probabilities for matching of the nodes of G 2 and G 3 

when the approximate Hermitian matrix is used. 



(c)Probabilities for matching of the nodes of Gi and G2 (d)Probabilities for matching of the nodes of Gi and G 2 

when the approximate Hermitian matrix is used. 

FIG. 2. Expected probability outcome for the alignment of pairs of the networks 


be obtained from both of the figures. 


VI. THE OVERALL ALGORITHMIC 
COMPLEXITY 

The computational complexity of a quantum algorithm 
is measured by the number of fundamental quantum 
gates-e.g. one and two qubit gates-in the circuit which 
simulates the algorithm. The complexity of PEA is dom¬ 
inated by the complexity of the circuit implementing the 
given operator A. Implementation of A order of N re¬ 
quires 0(N 2 ) number of quantum gates However, it 
is proven that the efficient (polynomial time in the num¬ 
ber of qubits) simulation of a sparse operator on quantum 
computers is possible when the number of entries in A is 
bounded polynomially in the number of qubits and the 
norm of the matrix is less than or equal to the degree 
of this polynomial [21:] . There have been also algorithms 
presented to simulate such sparse operators in polyno¬ 
mial time 0,113 • A s a result, since the adjacency ma¬ 


trices of the PPI networks sparse, they can be simulated 
efficiently on quantum computers. 

In the implementation of the operator A , quantum cir¬ 
cuits for each graph can be generated separately since it is 
the Kronecker product of the column-wise normalized ad¬ 
jacency matrices (normalized in the sense to make them 
stochastic): A = Ai ® A 2 0 ■ ■ ■ ® A m , where m represents 
the number of networks. Therefore, the total complexity 
of implementing A can be defined as: 

o(poly(log( |Vi|)) + poly(log(\V 2 \)) .. .poly(log{\V m \))'), 

( 8 ) 

or more concisely: 

o(m x poly(log(\V max \)) (9) 

where |I^n QX | defines the maximum number of nodes in a 
graph. This indicates an exponentially faster evaluation 
time for the implementation of A than the implementa¬ 
tions on classical computers. Note that if the adjacency 
matrices are not sparse, then the complexity for this part 
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becomes 0(rn|V maa; | 1 2 3 4 5 6 7 8 9 10 ) which still provides exponential 
speed-up in comparison to 0(|-E| m ) classical complexity. 
Moreover, in the case of the incorporation of additional 
data: i.e., using A+B instead of A ; the above complexity 
arguments hold by assuming A and B commute and B 
is efficiently simulatable. 


The complexity of the matching part of the algorithm 
is related to the number of measurements applied to the 
system. If only one kind of measurement setting (mea¬ 
sure a register alone and the rest together to draw the 
conditional probabilities) is used, then this part requires 
polynomial time (polynomial by the number of qubits) 
because it is related to the number of qubits. On the 
other hand, if one also uses different combinations of the 
registers in the measurements and consider them together 
to match the nodes (This can be only useful when addi¬ 
tional data is incorporated.), then the complexity to store 
the statistical results of the measurement outcomes and 
find matching from these outcomes may grow exponen¬ 
tially with the number of networks. However, as men¬ 
tioned in Sec lIVBI if no additional data is used, then the 
system is separable and all possible measurement settings 
produce the same output. Therefore, only one measure¬ 
ment setting can be used. 


VII. CONCLUSION 

In this paper, we have presented a quantum approach 
for the alignment of multiple networks by adapting quan¬ 
tum phase estimation algorithm. In particular, we have 
showed that the principal eigenvector of a stochastic ma¬ 
trix used in IsoRank algorithm for the alignment can be 
found exponentially more efficiently on quantum phase 
estimation algorithm. Since the final quantum state rep¬ 
resenting the eigenvector in the phase estimation algo¬ 
rithm is not classically available, adapting a conditional 
measurement scheme, we have also showed a matching 
algorithm to obtain the alignment result from this state. 
In addition, since the stochastic matrices are generally 
not Hermitian, we have also discussed how to approx¬ 
imate them for the simulation on quantum computers. 
Finally, we have used three simple networks and showed 
the numerical alignment results for them. While the ap¬ 
proach discussed here follows mainly IsoRank algorithm, 
we believe it shall also pave the way for the applications 
of other spectral alignment methods on quantum com¬ 
puters. 
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(a)Probabilities for matching of the nodes of G i, G 2, and G3. 
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FIG. 3. Probabilities for matching of the nodes Gi, G 2 , and G 3 . While the x-axis represents the nodes, aibjCk ; the y-axis is the 
unnormalized conditional probabilities: if \regd) is in a* state, then the probability to measure | regc 2 ) \ re 9G 3 ) in bjCk state 
in the collapsed quantum state is shown. 










